Saffman— Taylor instability in a non-Brownian suspension: 
finger selection and destabilization 
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We study the Saff man-Taylor instability in a non-Brownian suspension by injection of air. We find 
that flow structuration in the Hele-Shaw cell can be described by an effective viscosity depending on 
the volume fraction. When this viscosity is used to define the control parameter of the instability, 
the classical finger selection for Newtonian fluids is recovered. However, this picture breaks down 
when the cell thickness is decreased below approximatively 10 grain sizes. The discrete nature of 
the grains plays also a determinant role in the the early destabilization of the fingers observed. 
The grains produce a perturbation at the interface proportional to the grain size and can thus be 
considered as a "controlled noise". The finite amplitude instability mechanism proposed earlier by 
Bensimon et al. allows to link this perturbation to the actual values of the destabilization threshold. 

PACS numbers: 47.15.gp, 47.54.-r, 47.55.Kf, 83.80.Hj 
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Introduction - When a low viscosity fluid like air dis- 
places a viscous, immiscible fluid in a thin channel or 
Hele-Shaw cell, an instability develops at the interface, 
leading to the formation of finger like patterns, called vis- 
cous fingers or Saff man-Taylor instability [1]. Since the 
early work of Saffman and Taylor, this problem has re- 
ceived much attention not only because of its practical 
importance but also since it represents an archetype of 
many pattern forming systems [1, 0, 0] • Driven by prac- 
tical and fundamental interests, several viscous fingering 
studies have lately been extended to non-Newtonian flu- 
ids where a wide variety of strikingly different patterns 
are found [5]. 

Here, we use the Saff man-Taylor instability in a Hele- 
Shaw channel as a model system to study the dynamical 
properties of the interface between a pure fluid and a non- 
Brownian suspension. This kind of particle laden fluid is 
known to structure under shear as the particles migrate 
towards regions of lowest shear and thus change the local 
viscosity of the flow. As a consequence, this is a fluid 
with complex evolutive properties. A central question 
is to understand whether a simple effective continuum 
description remains valid and how pattern selection and 
finger stability can be affected by the fluid structuration 
as well as the discrete nature of the grains. 

Set-up and characterization of the suspensions - The 
experiments are performed in an Hele-Shaw cell of length 
1 m formed by two 1.5 cm thick glass plates separated by 
a thin mylar spacer (Fig. [1]). The cell thickness h can be 
varied (6=0.75 - 1.43 mm) as well as the width W of the 
channel {W=2 - 4 cm). The thin channel is initially filled 
with a non-Brownian suspension that is then displaced by 
air. The suspensions are formed by spherical polystyrene 
beads from Dynoseeds suspended in a Silicon oil. We use 
different grain diameters 1^=20, 40, 80 or 140 /im, their 
density is p=1050~1060 kg/m^. The grains are dispersed 
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FIG. 1: Schematic drawing of the experimental set-up. Evo- 
lution of a typical finger with increasing finger velocity (a)-(f). 



in a modified Silicon oil (Shin Etsu SE KF-6011) such as 
to obtain density matching at a value p=1070 kg/m^. We 
measured the viscosity of the pure fluid r^o = 191 mPa.s 
as well as its surface tension 7 =21±1 mN/m at 21 °C. 
The fact that the particles and the suspending fluid are 
density matched allows us to control the volume fraction 
over a wide range. Here we work with volume fractions 
below 40 % which allows to avoid jamming problems and 
wall slip occurring at higher volume fractions. 

A constant overpressure is applied at the inlet whereas 
the outlet is maintained at atmospheric pressure. The 
advancing finger tip is observed using a CCD camera 
mounted on a movable system which allows a manual 
tracking of its position. The camera is coupled to a mi- 
crocomputer for direct image acquisition. Note that in 
this configuration, the applied pressure gradient is not 
constant and therefore, the finger accelerates when prop- 
agating through the cell. We verified that this accelera- 
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FIG. 2: a) Relative viscosity of the suspensions (L)=80/im) as a function of the grain fraction obtained from rheological 
measurements riRh{(j)) /rj^ (□) and from flow in the Hele-Shaw cell ric{(l))/rio (•)■(-) theoretical prediction of Zarraga et al. [6|. 
b and c) Relative finger width A as a function of 1/B = 12{W/b)'^ {r]cU/-f) for the same suspensions. Cell width VK=4 cm and 
cell thickness 6= 1.43mm (b) and 6=0. 75mm (c) 



tion is slow enough and does not influence the observed 
finger properties (i.e. for a given tip velocity, we have no 
dependence on the applied overpressure). 

Finger width selection - For Newtonian fluids, the 
width w of the viscous fingers is determined by the capil- 
lary number Ca = r]U/j^ the ratio of viscous to capillary 
forces; U being the finger velocity, r] the viscosity and 7 
the surface tension of the fluid. More precisely, the rel- 
ative finger width X = w/W is di function of the dimen- 
sionless control parameter 1/B = 12{W/b)'^Ca^ which 
involves the cell aspect ratio. The mean flow (averaged 
over the cell thickness) is governed by Darcy's law which 
links the fluid velocity V to the pressure gradient and far 
away from the finger, reduces to: V = where 
VP is the applied pressure gradient. 

Therefore, by analogy with Newtonian fluids, we seek 
to establish the relevant control parameters of the injec- 
tion process. A first step will be to measure the suspen- 
sion viscosities as a function of volume fraction. Since 
we want to evidence the influence of structuration ef- 
fects due to the confinement in the Hele - Shaw cell, we 
will compare the viscosities obtained from the commer- 
cial rheometer to those extracted directly from Darcy's 
law in the Hele-Shaw cell. 

The suspension viscosities r] (for 0=40 and 80 fim) 
were obtained by rheological measurements using a dou- 
ble Couette geometry rheometer (Haake-RS600) of gap 
width 2*0.25 mm and mean radius 20 mm. The gap 
width being small with respect to the radius, the veloc- 
ity profiles in the gap can be considered as linear. In the 
range of shear rates tested (7=0.1 - 100 s~^) correspond- 
ing also to the typical shear rates of the Hele-Shaw exper- 
iments, the suspensions behave as a Newtonian fluid with 
a viscosity r]Rh{(t)) independent of the grain diameter D 
and well described by models used in recent literature (as 
for example Zarraga et al [6|]) (see Figure [2^ open sym- 
bols). These results confirm the absence of aggregation 
in the range of volume fractions considered. 

Now we investigate the rheology of our suspensions 
directly in the Hele-Shaw cell. To do so, we systemat- 



ically established a Darcy's law for all suspensions and 
cell geometries we have considered. The details of the 
procedure can be found in [3, 18|] • The results indicate the 
existence of an effective viscosity parameter specific for 
flow in a Hele-Shaw cell ?7c(0)- In all the studied con- 
figurations, this viscosity was found to be independent 
of the cell geometry (i.e. b and W) as well as the grain 
diameter D and is only function of volume fraction 0. 

We found for increasing cj) an increasing deviation of rjc 
from the corresponding rheometer viscosity rjRh (figure [2] 
a close symbols). The fact that r]c is lower than r]Rh can 
easily be explained by the effect of flow structuration due 
to confinement. It is well known that particles migrate to 
zones of low shear rate [9], here in the middle of the gap. 
Therefore, the flow profile should deviate from an ideal 
parabolic profile and evolve towards a flatter profile [10], 
independent of V and solely dependent on (j). A steady 
profil is only reached after a certain flow distance below 
which we do not consider our experiments. In addition 
we also verified that the grains in the suspension do not 
alter the surface tension. 

We systematically studied the selection of the finger 
width for stable fingers different values of the volume 
fraction (j). Typical results are displayed on figure [2)3 and 
c. The solid line represents the pure fluid measurements 
and thus corresponds to the classical result for Newtonian 
fluids. Importantly, we choose to use 7?c(</>) to define the 
control parameter 1/B = 12{W/b)'^{r]cU/j). We notice 
that this choice rescales our data when the ratio between 
the cell thickness and the grain diameter is larger then ap- 
proximatively 10. Below b/D ~ 10 significant deviations 
from the classical result appear: fingers are slightly larger 
for low 1/5, thinner for intermediate 1/B - this narrow- 
ing is more important for large volume fraction - and then 
eventually A will join the classical curve for larger values 
of 1/B. In all cases, for large 1/B the asymptotic value 
1/2 for the relative finger width seems to be recovered. 

So far, we were not able provide any fully satisfactory 
argument to explain the deviation from the classical re- 
sults. One might however speculate that effects like a 
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FIG. 3: Control parameter at the threshold 1 /Be as a function 
of the grain fraction cj). VK=4cm, D=80/im and 6=0.75 mm 
(□) and 6=1.43 mm (o). Inset: Localisation of the perturba- 
tion (o) during a typical experimental finger destabilization. 

depletion of grains in front of a receding meniscus 
(where 3D effects become important) would lead to a 
lower viscosity at the finger tip and thus induce a finger 
narrowing as observed for shear thinning fluids [12]. One 
could also consider other effects like normal stresses exist- 
ing in granular suspensio ns p] , that have been shown to 
increase the finger width [12'] and might thus eventually 
be responsable for the fact that the narrow fingers get 
wider again and stem to the classical results at a higher 
control parameter. 

Finger destabilization - When increasing the finger ve- 
locity and thus the control parameter 1/B one observes 
a destabilization of the fingers. First, the sides of the 
fingers start to ondulate and when the velocity is fur- 
ther increased, the finger gets unstable by tip splitting. 
A typical evolution of such a destabilization can be seen 
on figure [H Note that in some cases asymmetric fingers 
(not shown here) have been observed like predicted by 
Ben Amar et al. 

This destabilization scenario is also observed in clas- 
sical Saff man-Taylor experiments and has so far been 
attributed to the noise in the cell hence making it diffi- 
cult to predict exact values for this instability threshold 
[3, In the case of our suspension, we identified the 
instability threshold for different volume fractions and 
cell thicknesses. To have a reproducible and objective 
method of detection, we have measured by image pro- 
cessing the fluctuations of the finger width 5\ and detect 
when this quantity starts to increase. The onset of fluctu- 
ation growth defines the stability threshold and therefore, 
the critical control parameter 1/Bc. 

When the values oil/ Be are plotted against the grain 
fraction (see figure [3]) for two different set ups (different 
cell thickness, same grain size), we surprisingly find that 
as soon as a small amount of grains (as small as 1 %) 
is added, the value of the threshold drops strongly and 
thereafter varies weakly with the grain fraction (if not 



TABLE I: Critical control parameter 1/Bc for different ex- 
perimental configurations and cj) — 10%. 



N config. 


D 


h 


W 


l/B, 


1 


80 /im 


0.84 mm 


20 mm 


1210 


2 


20 /im 


0.75 mm 


40 mm 


3220 


3 


40 /im 


0.75 mm 


40 mm 


2670 


4 


80 /im 


0.75 mm 


40 mm 


1600 


5 


80 /im 


1.43 mm 


40 mm 


1530 


6 


140 /im 


1.43 mm 


40 mm 


1150 



constant within the experimental uncertainties). Note 
that we found different values for the threshold in pure 
fluid for the different configurations and thus two differ- 
ent intrinsic levels of noise. However in the presence of 
grains the value of the threshold is well defined. 

When looking closer at the destabilization scenario of 
a single finger (figure [3] inset), one observes that pertur- 
bations nucleate close to the finger tip and are then ad- 
vected to the side of the finger. A possible explanation for 
this destabilization mechanism was given by Bensimon et 
al. [4] . They propose that a perturbation would nucleate 
at the finger tip, where the normal velocity is the highest 
and where the growth rate (given by the linear stability 
analysis of the Saff man-Taylor problem [16]) is largest. 
While the finger continues to grow, the perturbation is 
advected to the side of the finger where the normal veloc- 
ity goes to zero and thus, the perturbation growth stops. 
In the process of advection, the perturbation is stretched 
and consequently, its amplitude is decreased. For a given 
control parameter l/^c, one needs a given, finite ampli- 
tude of the initial perturbation Ai to be able to obtain a 
perturbation with a given final amplitude Af when the 
side of the finger is reached. More quantitatively, Ben- 
simon et al. derived the following relation that describe 
the "finite amplitude instability": 

Af ^ Ai exp (0.106 1/5^^) . (1) 

In the following, we will test if the destabilization ob- 
served in our situation can be described by the previous 
mechanism and if the existence of grains in the viscous 
fluid can be directly linked to the instability onset. A first 
indication is that we observe a destabilization threshold 
hardly affected by an increasing grain fraction; the above 
described mechanism is indeed independent of the wave- 
length of the perturbation which one might be tempted to 
link to the grain fraction. In their theoretical approach, 
Bensimon et al 0] consider the final amplitude Af to be 
proportional to the channel width W . For our analysis, 
we directly use the value of the fluctuations of the finger 
width observed at 1/Bc. Those are of course via the 
finger width proportional to W and we consider a sinu- 
soidal form of the fluctuations leading to Af = 8\W2\/2. 
A natural assumption is to relate the amplitude of the 
initial perturbation Ai to the grain diameter D\ Ai (x. D. 
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FIG. 4: Experimental values (□) and theoretical prediction 
form Bensimon et ai (-). The numbers correspond to the 
experimental configuration given in table ID Additional ex- 
periments (•) using a different silicon oil (Dow Corning 704). 

In table [J we sumarize the results obtained for different 
grain sizes D and different cell widths W. When we plot 
^JXjBc as a function of ln(A//I)) obtained for ^ = 10%, 
different cell geometry and different grain sizes (Fig. [4]), 
we obtain a very satisfying quantitative agreement with 
the theoretical result of Bensimon et al when the initial 
perturbation is taken to be: Ai ^ D/20. 

Now, to provide a rational justification for our choice 
for the initial perturbation, we need to show that a parti- 
cle approaching a free interface at a stagnation point flow, 
as is the case of the finger tip, would indeed be able to de- 
form this interface with the correct amplitude. Hoffman 
and Montiel et al. [18] established a link between the 
free interface amplitude of perturbation and a capillary 
number defined on the scale of a particle: Cap = ^^Ca. 
In our case, the range of capillary numbers would be 
(7ap=5*10~^ - 2*10~^ and thus, the corresponding val- 
ues for the amplitude of the perturbation are of order of 
0.02 — 0.1 X D. Therefore, this is indeed in good agree- 
ment with our choice of the amplitude of the initial per- 
turbation and it confirms the validity of the theoretical 
prediction of Bensimon et al [4] . This is a striking emer- 
gence of the particulate nature of the suspension, the 
grains acting here as a controllable noise amplitude. 

Conclusion - We have shown that the characteristics 
of the stable Saffman - Taylor fingers invading a non- 
Brownian suspension, are mainly determined by the ef- 
fective properties of the suspension in the cell far from 
the fingers. The structuration effects of the confined par- 
ticulate flow can simply be taken into account by an ef- 
fective viscosity, which allows to rescale the width selec- 
tion curve onto the classical results. However, when the 
grain size becomes important compared to the cell thick- 
ness, systematic deviations from the classical result are 
evidenced. 

We also studied the destabilization of the Saffman- 
Taylor fingers and we showed that individual grains per- 



turbe the interface between the air and the suspension 
and lead to a premature destabilization of the fingers. 
The threshold of instability is found to match quantita- 
tively the theoretical prediction of Bensimon et al. To 
our knowledge this is the first time that an experiment 
allows to control the initial "noise" in the cell and thus 
to investigate a mechanism of finite amplitude instability 
directly. 

Recently a number of studies have reported oscillations 
of the finger width, observed for example for low capillary 
number [19] or fixed perturbations of the cell thickness 
(2q| . A more close characterization of the oscillations re- 
sulting from the premature destabilization in our system 
might reveal similarities between the different systems. 

We wish to thank Michel Cloitre and Fabrice Monti 
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